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I .  INTRODUCTION 


During  the  past  15  vears  a  number  of  powerful  high-frequency  fHF) 
radio  transmitters  have  been  developed  solelv  to  transmit  waves  strong 
enough  to  produce  measurable  nonlinear  effects  in  the  ionosphere.  At 
first  the  onlv  goal  was  to  produce  a  detectable  change  in  electron 
temperature,  so  the  early  tests  were  called  heating  experiments.  Manv 
-vpes  of  nonlinear  effects  were  observed,  however,  and  such  tests  are 
now  called  ionospheric  modification  experiments.  Several  authors  ( 1- 
3 |  have  summarized  the  phenomena  observed  in  modification  experiments. 

Although  some  Soviet  work  [4-7]  has  used  oblique  incidence,  vir¬ 
tually  all  modification  experiments  have  transmitted  the  powerful 
modifying  wave  at  vertical  incidence.  Radar  and  communication  systems 
operate  at  oblique  incidence,  however.  For  reasons  given  below,  an 
obliquely  incident  modifying  wave  might  not  produce  the  same  effects 
as  a  vertically  incident  one.  It  is  therefore  of  interest  to  examine 
the  effects  of  a  strong  oblique  wave  on  the  ionosphere  and  on  waves 
that  propagate  through  the  modified  region. 

The  nonlinear  processes  that  lead  to  ionospheric  modification  are 
divided  into  two  categories:  (1)  simple  heating,  which  alters  tem¬ 
perature  and  hence,  reaction  rates,  collision  frequencies,  and  par¬ 
ticle  densities,  and  (2)  generation  of  parametric  instabilities,  which 
cause  a  myriad  of  phenomena,  including  spread-F  and  geomagnetic-field- 
aligned  irregularities.  Both  categories  can  change  the  amplitude, 
phase,  or  path  of  a  signal  that  traverses  the  modified  region,  and 
both  depend  strongly  on  the  incidence  angle  of  the  modifying  wave. 

As  the  incidence  angle  is  increased,  heating,  which  is  propor¬ 
tional  to  the  local  power  density,  is  affected  by  two  competing 
trends.  The  field  is  weakened  by  the  increased  path  length,  but 
strengthened  by  focusing  near  caustics  [8].  Field  and  Warber  [9,  10] 
showed  that  focusing  near  caustics  overcomes  the  greater  geometric 
spreading  suffered  bv  oblique  modifying  waves,  and  that  intense 
electric  fields  can  occur.  They  only  calculated  electric  field 
strength,  however,  and  did  not  account  for  transport  processes;  nor 
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did  they  calculate  the  electron  density  changes  caused  bv  chemistry 
changes  and  diffusion  driven  by  such  focused  fields.  Most  parametric 
instabilities  cannot  be  excited  by  oblique  waves,  because  only  verti¬ 
cal  waves  can  reach  heights  where  matching  conditions  on  frequency  are 
me  t . 

The  present  report  extends  the  work  of  Field  and  Warber  to  include 
the  effects  of  heat  conduction  and  diffusion.  Those  processes  are 
more  important  for  oblique  incidence  than  for  vertical  incidence  [11], 
because  at  oblique  incidence  the  intense  fields  are  strongly  localized 
near  caustics,  which  are  narrow.  The  resulting  gradients  are  stronger 
than  those  caused  bv  vertical  waves,  whose  energy  is  deposited  over 
wider  ionospheric  regions.  This  report  emphasizes  the  underlying 
theory  and  derives  equations  suitable  for  numerical  calculation.  In 
addition,  sample  cases  are  presented  that  illustrate  how  heat  conduc¬ 
tion  and  ambipolar  diffusion  affect  temperature  and  density  changes  in 
the  caustics  of  intense  oblique  waves. 
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II.  EQUATIONS  GOVERNING  ELECTRON  HEATING  AND  TRANSPORT 


The  equations  that  govern  the  changes  in  electron  temperature  and 
density  produced  by  an  intense  radio  wave  in  the  ionosphere  are 
presented  in  this  section.  Analytic  solutions  for  idealized  situa¬ 
tions  that  are  useful  for  understanding  the  physical  processes  are 
provided . 

Some  of  the  equations  below  are  similar  to  those  in  Gurevich  [1], 
whose  notation  we  use  whenever  convenient.  Certain  lengthy  deriva¬ 
tions,  not  given  bv  Gurevich,  are  outlined  in  the  appendixes. 

Because  the  magnetic  pressure  greatly  exceeds  the  plasma  pressure 
at  altitudes  of  interest,  we  consider  only  field-aligned  particle 
motion  and  energy  transport,  ignoring  the  small  cross-field  drift  and 
thermal  conduction.  We  use  a  coordinate  system  with  the  x-direction 
along  the  geomagnetic  field.  That  coordinate  is  related  to  height  z 
bv  x  sin  0  =  z,  where  0  is  the  dip  angle. 


ELECTRON  TEMPERATURE 

If  we  neglect  ion  heating,  the  electron  temperature  T  is  described 
bv  the  following  equation  (e.g.,  Meltz  and  LeLevier  [12]): 
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where  T0  is  the  ambient  electron  temperature.  E  is  the  electric  field 
of  the  modifier  wave,  S  is  the  fractional  energy  lost  by  an  electron 
per  collision,  and  i/e  the  electron  collision  frequency. 

The  quantity  Ep  in  Eq .  (1)  is  defined  as 
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where  K  is  Boltzmann's  constant,  m  is  the  electron  mass,  e  is  the 
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electron  charge,  and  ui  is  the  modifier  wave  angular  frequency.  Be¬ 
cause  Ep  has  units  of  volts  per  meter,  it  is  often  called  the  charac¬ 
teristic  plasma  field.  The  form  given  by  Eq .  (2)  is  valid  when  the 
wave  frequency  greatly  exceeds  the  electron  gvrof requencv . 

The  quantity  k  in  Eq .  (1)  is  the  thermal  diffusivity  and  is  given 


k  ~  A 


where  A  is  a  constant  that  depends  on  the  sophistication  of  the  theory 
used  to  derive  k.  .  Although  A  can  be  as  large  as  2.  it  is  often  as¬ 
sumed  to  be  unitv.  The  thermal  diffusivity  accounts  for  heat  conduc- 
t  ion . 

The  coefficients  k  and  ue  both  depend  on  temperature,  so  Eq.  (1; 
is  non  linear.  Moreover,  the  dependence  of  k.  on  x  prevents  Eq .  il) 
from  assuming  the  form  of  a  canonical  diffusion  equation.  The 
criterion  for  linearizing  Eq .  (1)  is  that  AT  =  T  -  T0  «  T .  which  is 
always  true  if  E2/E2  is  small,  and,  in  the  presence  of  strong  heat 
conduction  can  be  true  even  if  E2/E2  is  of  order  unity.  The  x- 
dependence  of  the  thermal  diffusivity  k  can  then  be  ignored  if  the 
gradients  in  the  ambient  ionosphere  are  weaker  than  those  caused  by 
the  modifying  wave.  We  make  both  of  those  simplifications  in  the 
present  analysis,  and  rewrite  Eq .  (1)  in  the  form 


2  [21 
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dx  E  o 

p 


where,  following  Gurevich  [1],  we  have  defined 


rT  =  (&l>  )  1  s 
e 


. 2  2 

L.p  =  m 


and  y  is  computed  using  T  =  TQ .  The  quantity  r  is  the  characteristic 
time  required  to  achieve  most  of  the  temperature  change.  The  quantity 


Lp  is  the  characteristic  distance  over  which  heat  is  conducted  in  time 
r ^ ,  and  thus  measures  the  dimension  of  the  region  over  which  the 
energy  deposited  bv  the  modifying  wave  is  spread  through  heat  conduc- 
t  ion . 

If  we  use  locaL  height  averages  for  rp  and  Lp ,  which  is  equivalent 
to  assuming  that  those  coefficients  are  locally  he ight- independent . 

Eq .  (4)  can  be  solved  as  shown  in  Appendix  A.  The  result  is: 
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Once  the  field  E  of  the  modifying  wave  has  been  calculated.  Eq .  (7) 
can  be  used  to  calculate  both  the  time  and  space  dependence  of 
electron  temperature. 

The  steady-state  reached  by  the  electron  temperature  after  the 
initial  transient  period  may  be  of  interest.  Appendix  B  shows  that 
the  asymptotic  behavior  of  Eq .  (7),  applicable  to  late  times  where 
t  »  rp.  is  simply 
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ELECTRON  DENSITY 

If  we  again  use  Locally  height-averaged  values  for  ambient  iono¬ 
spheric  properties,  we  can  write  the  equation  that  governs  the 
electron  density  in  the  following  form: 
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In  Eq .  (9),  P  is  the  ion-pair  production  rate,  a  is  the  eLectron-ion 
recombination  coefficient,  and  /3  is  a  loss  rate  that  accounts  for 
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The 


interactions  between  the  ion  0+  and  the  molecules  N2  and  (>2  . 
ambipolar  diffusion  coefficient  is 


D 
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and  the  ambipolar  thermal  diffusivity  is 


KT 
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where  R^  =  T0/(T0  s-  T^),  M  is  the  mean  ion  mass.  Tj^  is  the  ion  tem¬ 
perature,  and  i/£  is  the  ion-neutral  collision  frequency.  Equation  (9) 
thus  accounts  for  diffusion  caused  by  excess  particle  concentrations, 
diffusion  driven  by  thermal  gradients,  and  density  changes  brought 
about  through  the  temperature  dependence  of  reaction  rates. 

The  fractional  density  change  is  always  small,  so  we  can  writ** 

N  =  N  +  n  ,  (12) 

o 


where  N0  is  the  ambient  electron  density.  We  can  also  use  the  fact 
that  n  «  N0  to  linearize  Eq .  (9).  The  result  is: 
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where  a0  a  a(T0)  is  the  unperturbed  recombination  coefficient  and, 
because  the  field  E  is  not  strong  enough  to  alter  the  ambient  produc¬ 
tion  rate  P,  we  have  used 


P  =  a  N2  +  0  N  .  (14) 

00  00 

Because  fl  is  virtually  independent  of  temperature,  we  have  neglected  a 
term  proportional  to  0  -  fl0.  Bv  defining 
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(15) 
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-■e  can  rewrite  Eq .  (13)  in  a  form  analogous  to  Eq .  (4) 
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We  reconcile  Eq .  (17)  with  the  definition  of  y  used  by  Gurevich  [1)  by 
noting  that  the  electron  loss  rate  qr  is  qr  =  aN2  +  /3N ,  so  = 
(3qr/5N)-^  and  da/8T  =  N~2aqr/aT.  The  time-dependent  solution  to 
Eq .  (18)  car.  be  written  in  terms  of  a  Green's  function,  but  that 
result  is  too  complicated  to  be  useful.  Therefore,  we  restrict  atten¬ 
tion  to  the  following  steady-state  solution,  which  is  derived  in 
Appendix  C' 
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As  we  have  shown  for  the  temperature  change  AT,  once  the  modifying 
field  E  is  calculated  using  a  method  described  by  Field  and  Warber 
[9],  the  density  change  n  can  be  calculated  by  performing  a 
straightforward  numerical  integration. 
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III.  VALUES  OF  TRANSPORT  COEFFICIENTS  AND  ELECTRON  LOSS  RATES 


The  coefficients  and  reaction  rates  that  appear  in  Eqs.  (4) 
through  (8),  (18),  and  (19)  are  presented  in  this  section.  Most  of 

those  quantities  cannot  be  specified  precisely,  because  either  they 
depend  on  ionospheric  properties,  which  fluctuate,  or  thev  depend  on 
processes  that  are  poorly  understood.  Gurevich  [1]  suggests  nominal 
values;  we  consider  them  reasonable  and  use  them  as  inputs  in  our 
analysis.  The  following  discussion  gives  the  rationale  for  using 
those  values,  and  indicates  the  uncertainties.  We  limit  discussion  to 
altitudes  between  200  km  and  300  km,  because  that  is  where  the  most 
important  heating  occurs. 


MODEL  IONOSPHERE 

Table  1  summarizes  the  ionospheric  parameters  suggested  bv 
Gurevich  (1).  Those  parameters  represent  nominal  values  suitable  for 
examples,  but  their  values  vary  considerably,  depending  on  sunspot 
activity,  latitude,  and  local  time.  In  our  notation,  N^  denotes  the 
number  density  of  nitrogen  molecules,  Nnq+  denotes  the  density  of 
nitric  oxide  ions,  and  so  forth.  Oxygen  atoms  are  the  main  neutral 
constituent  at  the  altitudes  shown. 


RESPONSE  TIMES,  TRANSPORT  COEFFICIENTS,  AND  CHARACTERISTIC  LENGTHS 
Nominal  values  for  response  times  rj  and  r^,  transport  coef¬ 
ficients  k.  and  Da,  and  characteristic  lengths  Lj  and  LN  are  derived  by 
inserting  parameter  values  from  Table  1  into  formulas  given  in 
Sec.  II.'*  The  resulting  values  are  listed  in  Table  2.  Although  the 


"No  simple  formula  is  given  for  the  energy  loss  fraction  A.  be¬ 
cause  different  processes  contribute  to  6  at  different  altitudes.  At 
altitudes  so  high  (>  250  km)  that  elastic  collisions  between 
electrons  and  ions  predominate,  A  is  simply  3m  H.  At  lower  altitudes, 
where  inelastic  collisions  with  molecules  are  important.  A  must  be 
calculated  from  the  appropriate  energy-dependent  cross  sections. 
Gurevich  (1|  gives  altitude-dependent  recipes  for  calculating  A. 
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TABLE  1.  Model  Ionosphere 


10 


SOURCE. 


TABLE  2.  Transport  and  Loss-Rate  Coefficients. 


250  25.0  *00  370  0.61  103  16  0.52 


,  i  -  .  A  ia /-it  c:  c  *  '  :■»  /1i  f  fnronf  On^*^  hf* 

nominal  values  given  au  userui  toi  cuscus^w.. 

defined  Co  correspond  Co  the  actual  model  ionosphere  used  for  a 
specific  calculation. 

ELECTRON  LOSS  RATES 

Gurevich  (1)  suggests  the  following  formulas  for  the  loss  rates  a 
and  0: 


V*  °3 

•  "i  ^  2  n 


whe  re 


-  5  x  10 


-7  300 


(21a) 


„2  -  2.2  X  10-;  m 


(21b) 


0  =  10'12  N  +  2  X  10"U  N0 


<210 


Other  authors  give  somewhat  different  values  for  the  above  loss 
rates.  For  example.  McEwan  and  Phillips  [13]  indicate  that  aj  behaves 
as  (300/T)^-°.  whereas  Bates  [14]  uses  the  relationships 


=  4.2  x  10 


-7  300  1 
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L . 6  x  10 


300 

T 


(23) 


where  0.37  <  <  0.9  and  0.55  <  Vo  £.  0.7.  We  see.  therefore,  that 

Gurevich  uses  a  somewhat  stronger  temperature  dependence  tor  recom¬ 
bination  than  suggested  by  other  authors,  particularly  for  recombina¬ 
tion  with  nitric  oxide.  The  actual  dependences  are  uncertain,  but. 
Gurevich's  formula--which  we  use--predicts  a  density  perturbation  n 
that  is  on  the  high  side  of  the  uncertainty  range. 
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IV.  SOLUTIONS  FOR  SOME  SPECIAL  CASES 


To  illustrate  the  dependence  of  electron  temperature  and  density 
on  the  various  coefficients,  solutions  are  presented  for  idealised 
cases.  The  equations  of  Sec.  II  can  be  integrated  for  these  cases 
analytically . 

LIMITING  CASE:  NO  HEAT  CONDUCTION  OR  DIFFUSION 

If  heat  conduction  and  diffusion  are  neglected,  which  is  valid  if 
the  dimension  of  the  heated  region  far  exceeds  the  lengths  Lpj  and  bp, 
then  L-p  =  L^  «  0  and  Eqs.  (4)  and  (18)  become 


•(2a  N  +  0)n  -  [a(T) 
o  o 


a(T  )]N 
o  o 


The  solution  to  Eq .  (24), 


exp ( - t/r  ) 


which,  for  late  times,  becomes 


r  t  »  t  2 

o  T  E 

P 


Thus,  in  the  absence  of  diffusion  or  heat  conduction,  the  relative 
temperature  change  is  simplv  the  squared  ratio  of  the  modifying  field 
E  to  the  characteristic  field  Ep .  In  that  limit,  the  heated  region  is 
"frozen"  and  has  spatial  features  identical  to  those  of  E‘ . 


In  the  same  steady-state  limit,  the  solution  to  Eq .  (25)  is 


N^[a(To)  -  a (T) ] 

t  »  r  2a  N  +  0 

N  oo 


(28) 


The  important  result  of  Eq .  (28)  is  that  if  there  were  no  diffusion, 
the  electron  density  perturbation  n  would  occur  solely  through  the 
temperature  dependence  of  the  recombination  coefficient  a.  Because 
recombination  is  unimportant  above  about  200  km,  virtually  the  entire 
electron  density  perturbation  at  such  high  altitudes  depends  upon 
ambipolar  diffusion,  omitted  from  Eq .  (28).  Stated  differently,  at 
altitudes  much  above  200  km  there  would  be  no  perturbation  if  there 
were  no  diffusion. 


SPATIALLY  UNIFORM  HEATING  PULSE 

If  the  heating  pulse  has  a  "square-wave"  shape  given  by, 


' 

E  |  x  |  s  a 

E(x)  =  •  °  (29) 

0  |x|  >  a  , 

.. 

where  E0  is  constant  and  2a  is  the  width  of  the  heated  region,  then 
Eqs.  (7),  (8),  and  (19)  can  be  integrated  straightforwardly. 

It  is  convenient  to  normalize  the  fractional  temperature  and 
density  changes  to  their  values  in  the  absence  of  transport  processes 
and  recombination.  We  therefore  write 

AT  =  (AT/To)/(E2/Ep)  ,  (30) 

and 

n  =  (n/NQ)/(E2/E2)  .  (31) 
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Without  heat  conduction,  diffusion,  or  a  temperature-dependent  recom¬ 
bination  rate,  the  normalized  temperature  and  density  perturbations 
are 


AT  =  1  , 
e 


n  =  0 


x  <  a  , 


(32) 


AT  =n=0  x  >  a 
e  1  1 


If  we  use  the  above  form,  Eq.  (7),  for  the  time-dependent  tempera¬ 
ture  ,  becomes 


-t/r. 


AT  =  c  +  - 
e  o  2 


-b+/t 

'  So  erf  1  f  +  So  I  [w(iC+)  +  w(ic; 


)] 


-  erf 


Vi'! 


■b  /t 


[w(iC_)  +  w(if^) ] 


(33) 


where  w  is  the  error  function  for  complex  arguments  [15)  and 

12 


b,  =  f_ 
±  T 


x  ±  a 


2LT 


r±  =  7 t/rT  +  yb±/t  , 


r±  =  -  y t/rt  +  yb±/t 


c  =  < 

o 


3 

O 


1  |x|  <  a  , 

.  0  I  x  |  >  a  , 

-1  | x i  s  a 

1  I  x  I  >  a 


16 


Equation  (8),  for  the  stationary  perturbation  temperature,  becomes 


ATg  =  G(x ,  Lt) 


(34) 


where 


' 

> 

(“'  hj  *  %  +  1  s„ 

..  |x|  -  a 

-  exp 

x  |  +  a 

”  h  J 

lt 

or 


■a/LT 


cosh  x/L_ 


<  a 


G(x,  Lt) 


-  I x I /lt 

e  sinh  a/L^,  |x|  >  a  . 


Equation  (19),  for  the  stationary  electron  concentration,  becomes 


+ 


RjP  +  7 


G(x,  Lx)  -  G(x.  Ln)' 

__ 


(33) 


where 


2  2 
'  -  VLx 


Equation  (35)  is  undefined  for  p  =  1.  In  this  case,  for  |x|  <  a, 


n  =  7G 


(x.  Ln]  +  (V  ♦  7) 


"a/LN 

~2 


7—  sinh  -  r—  cosh 
•  N  LN  LN 


(36) 


and  for  |x|  >  a, 
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(37) 


n  =  7g(x,  I^j  +  |rtp  +  7] 


-  M/i* 


|x|  a  a  ,  a 

d — L  sinh  - —  -  - —  cosh  - — 

Si  Si  Si  Si 


Peak  Temperature  and  Density  Perturbation  in  Limits  of  Weak  and  Strong 
Transport 

The  peak  values  of  AT  and  n  usually  occur  at  the  origin  (x  =  0)  . 
Setting  x  =  0  in  Eqs .  (34)  and  (35)  and  reverting  to  temperature  and 
density  that  are  not  normalized,  we  find 


and 


T  x  =  0  _2 
o  E  I 

P 


(38) 


N  x  =  0  _2 
o  E 

P 


RT^  + 


7Lt 


lt 


(e'a/L» 


(39) 


Equation  (38) 
by  the  factor 


shows  that  heat  conduction  reduces  the  peak  temperature 
(.  - 


The  limit  of  weak  heat  conduction  and  diffusion  is  reached  when 
a  »  Lfj  and  Lf,  which  can  occur  either  if  the  heating  pulse  is  very 
wide  (large  a)  or  the  diffusion  and  conduction  lengths  are  small 
(small  L) .  In  that  limit,  Eqs.  (38)  and  (39)  become 


and 


AT  _ ( 

T  a  »  L 
o 


n 


N 

o 


a  »  L 


(40) 


(41) 
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which  agree  with  Eqs .  (27)  and  (28).  From  Table  2,  we  know  that  7  is 
negligible  above  200  km. 

Throughout  much  of  the  F-region,  air-chemistrv  loss  rates  are  so 
small  and  insensitive  to  temperature  that  7  can  be  assumed  to  equal 
zero,  in  which  case  Eq .  (39)  becomes 

n 
N 

o 


rtln 


1/ln 


-a/L_ 


-  e 


(42) 


To  examine  the  limit  of  strong  transport,  we  assume  L[g  and 
L-p  »  a,  in  which  case  Eqs.  (38)  and  (39)  become 


n 

aRT 

N 

0 

lt 

LT  +  LN. 

(43) 


(44) 


=  R„ 


LT  +  ln 


AT 

T 


(45) 


where,  from  Table  2,  Ljq  <  Lj.  Equations  (44)  and  (45)  show  that 
strong  heat  conduction  reduces  the  peak  temperature  perturbation  by 
the  factor  a/Lj.  Moreover,  n/N0  is  small  whether  transport  processes 
are  weak  or  strong,  because: 

(a)  If  transport  processes  are  weak,  then,  according  to  Eq .  (41), 
(n/No)/(AT/TQ)  =  7  «  1. 

( b )  If  diffusion  is  strong,  then  n/N^  <  AT/Tq .  per  Eq .  (45). 

But,  according  to  Eq.  (43).  AT/T^  itself  is  then  of  order 
a/L^,,  which  is  small. 
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Transport,  which  is  necessary  to  create  a  density  perturbation  in  much 
of  the  F-layer,  therefore  suppresses  the  temperature  increase  AT  that 
drives  the  density  perturbation.  Although  large  values  of  AT  can  be 
created  if  the  heating  width  2a  is  great  enough,  the  fractional  den¬ 
sity  change  n/N0  is  always  much  smaller  than  unity. 

Numerical  Examples 

To  illustrate  the  temperature  behavior,  we  choose  Lp  =  5  km,  a 
somewhat  lower  value  of  Lp  than  occurs  for  the  F-region.  Figure  1 
shows  the  temperature  perturbation  caused  by  uniform  heating  for  10  s 
in  a  region  40  km  wide  (a  =  20  km).  The  thermal  response  time  is  r = 
5  s.  After  2  s,  the  interior  of  the  heated  region  has  a  nearlv 
uniform  temperature,  because  it  is  wider  than  the  conductivitv  length 
Lp,  and  has  not  been  heated  long  enough  to  produce  gradients  large 
enough  for  heat  conductivity  to  alter  the  temperature  appreciably. 
Later,  as  temperature  gradients  increase,  the  high  temperatures  at  the 
edges  of  the  heated  region  are  smoothed  by  conduction.  The  total 
heating  time  is  about  twice  the  thermal  response  time,  and  the  central 
temperature  attains  85  percent  of  its  stationary  value  after  10  s. 

The  maximum  stationary  value  of  the  temperature  perturbation  is  about 
the  same  as  would  be  achieved  without  conductivity,  because  the  con¬ 
ductivity  length  is  smaller  than  the  width  of  the  heated  region. 

Figure  2  illustrates  the  effect  of  beam  width  on  electron  temper¬ 
ature.  In  this  case,  Lp  =  17  km,  which  is  the  actual  conductivity 
length  for  200-km  altitude  during  the  day  (Table  2).  Figures  2a-c . 
show  the  temperature  change  resulting  from  uniformly  heated  regimes 
that  are  20- ,  5-,  and  1-km  wide.  The  final  temperature  perturbation 
is  much  smaller  than  that  attained  in  the  nonconductive  limit  because 
even  the  widest  beam  is  comparable  in  dimension  to  the  conductivity 
length.  Thus,  for  all  cases  shown,  the  heated  region  is  too  narrow 
relative  to  Lp  to  sustain  the  full  temperature  perturbation:  the 
narrower  the  heater  beams,  the  more  the  temperature  per turba t ion 
decreases.  For  the  smallest  (1  km)  heating  width,  the  peak  tempera¬ 
ture  perturbation  is  less  than  5  percent  of  that  in  the  nonconducting 
limit.  Table  2  shows  that  Lp  during  the  day  at  200-km  altitude  is 
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less  chan  chat  at  any  other  F-region  altitude,  day  or  night.  The 
normalized  temperature  perturbation  therefore  would  be  even  smaller  at 
greater  altitudes  or  at  night,  than  at  200  km  during  the  day. 

Figure  3  shows  the  electron  temperature  perturbation  at  low, 
medium,  and  high  F-region  altitudes  for  a  20-km  wide  heating  beam 
during  the  day.  Parameter  values  are  given  in  Table  2.  As  previously 
stated,  the  perturbation  is  greatest  at  200  km.  As  the  conductivity 
length  increases  with  altitude,  the  amplitude  of  the  temperature 
perturbation  decreases.  Table  2  shows  that  the  thermal  response  time 
increases  with  altitude,  because  collision  frequency  decreases  with 
altitude.  It  thus  takes  the  temperature  longer  to  reach  its  station¬ 
ary  value  as  altitude  increases. 

The  main  question  is  whether  a  heater  can  produce  electron  density 
gradients  that  are  large  enough  to  focus  or  defocus  an  incident  beam 
and  alter  the  received  sky  wave  at  a  ground-based  receiver.  Figure  4 
shows  the  daytime  stationary  temperature  and  density  perturbations  at 
200-km  altitude  for  three  values  of  width  2a.  At  200-km,  an  increase 
in  temperature  causes  a  decrease  in  the  recombination  rate  which,  in 
turn,  causes  an  increase  in  electron  density  n. 

Two  diffusion  processes  determine  the  electron  density  change: 

(1)  thermal  diffusion,  where  electrons  of  high  thermal  velocity  leave 
the  heated  region  (equivalent  to  thermal  expansion):  and  (2)  ambipolar 
diffusion,  where  electrons  in  regions  of  high  density  diffuse  to 
regions  of  lower  concentration.  The  two  processes  oppose  each  other 
--that  is  one  reason  why  density  changes  are  small. 

Figures  5  and  6  show  the  stationary  daytime  temperature  and  den¬ 
sity  changes  at  250  km  and  300  km,  respectively,  during  the  day  for 
1-,  5-,  and  20-km  heated  widths.  The  recombination  rate  is  little 
affected  by  temperature  changes,  so  y  ~  o  above  200  km  and  high  tem¬ 
peratures  produce  a  density  reduction  due  to  thermal  diffusion. 
Accompanying  the  concentration  deficit  within  the  heating  region  is  a 
concentration  excess  in  the  surrounding  region.  The  excess  concentra¬ 
tion  is  required  for  overall  conservation  of  electrons  because,  in  the 
absence  of  ion-pair  production  or  recombination  changes,  there  can  be 
no  net  change  in  the  number  of  electrons. 
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Stationary  temperature  (solid)  and  electron  density  (dashed 
changes  at  200-ka  altitude  during  day  for  three  heated  widths 
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V.  DISCUSSION 


The  transport  and  chemical  equations  that  describe  electron  den¬ 
sity  and  temperature  changes  produced  by  electron  fields  in  the  iono¬ 
sphere  have  been  derived  and  cast  into  a  form  suitable  for  numerical 
integration.  By  using  those  equations  in  conjunction  with  previously 
developed  methods  for  calculating  electric  fields  near  caustics,  the 
density  and  temperature  changes  produced  by  an  oblique  ionospheric  HF 
heater  can  be  computed  straightforwardly.  The  necessary  inputs  are: 
(1)  heater  power,  (2)  heater  radiation  pattern,  (3)  heater  frequency, 
and  (4)  model  ionosphere. 

Analytic  solutions  for  idealized  heater  beam-shapes  indicate  that 
heat  conduction  and  diffusion  are  the  dominant  processes  above 
220  km.  Those  processes  suppress  the  temperature  substantially, 
because  the  caustics  are  narrower  than  conduction  lengths,  but  are 
necessary  to  produce  density  changes,  which  influence  HF  propagation. 
The  fractional  density  changes  are  always  small.  At  altitudes  of 
220  km  and  below,  the  dominant  process  is  the  heating- induced  reduc¬ 
tion  of  the  electron-ion  recombination  rate  which  causes  an  increase 
in  the  electron  density.  For  that  reason,  greater  density  changes  can 
be  produced  at  200  km  than  at  higher  altitudes. 
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Appendix  A 
SOLUTION  TO  EQ.  (4) 


The  equation 


3AT  2  <32AT 

rT  at  LT  ,  2 
3x 


which  describes  the  temperature  perturbation  of  a  plasma  due  to 
the  electric  field  E,  is  of  the  general  form: 


equation:  u^  -  ku^  +  -yu  =  Q(x,t)  , 

boundary  conditions:  u(x,o)  =  f(x)  and  u  -*•  o  as  x  -*  ±®  , 


where 


u 


t 


u 

XX 


3u 

at 


We  simplify  the  equation  by  multiplying  it  and  its  boundary 
condition  by  eTc,  and  let  u  =  e“Ttw(x,t).  Then, 


so , 


and. 


-*yt 

H  -  ku  +  yu  =  e  (w  -  yw  -  kw  +  yw) 
t  XX  c  XX 


Q(x.t) 


y  t 

w-kw  =e  Q(x,t) 
t  xx 


R(x, t) 


w(x,o)  =  u(x,o)  =  f (x) . 


We  claim  that  the  solution  to  Eq .  (A.l)  is  as  follows: 


(A.l) 
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ft 


w(x.t)  = 


ds 


d£ 


{ (x-£ )  /4k(  t  -  s)] 
74 irk  (t  -  s) 


R(C.s)  + 


■[(x-O  /4kt] 


d£ 


74irkt 


(A. 2) 


where  R(£,s)  =  e7S  Q(£,s); 


hence , 


— y  £ 

u(x , t)  =  e  w(x . t) 


ds 


0 


t  r*  .  .  (x 

exp  -7<t  -  s)  -  — 

d*  - L - - zjj.  Q(jeiS) 


-<t  -  j)J 


7  4nk(t  -  s) 


+  e 


-7t 


d£ 


-[(x-O  /4kt] 


74irkt 


£(0  ■ 


PROOF  OF  EQ.  (A. 2) 
Let 


w-kw  =  R(x,t)  , 
t  xx 

w(x,o)  =  f(x)  , 

3w 


w 

XX 


We  assume  that  w,  R,  f,  and  their  second  derivatives  are  square- 
integrable  on  (-<»,  <»)  .  Let 


w(A , t) 


f°° 

-iAx  .  .  . 

e  w(x .  t )  clx 

*  -00 


f(0 
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denote  the  Fourier  transform  of  w(x,t).  Then, 


,  >  v  -lAx  ,  ,  , 

w  (  A  ,  t )  =  e  w  (x  .  t )  dx 
t  t 


f°° 

e  £kw^  +  R(x.t)j 

*  — nr> 


■  -lAx  ,  .  -lAx  „  , 

k  •  w  e  dx+  e  R(x.t)  dx  . 

xx 


Assuming  Lim  w(x.t)  =  0,  and  lim  wx(x,t)  =  0,  the  first  term  can  be 
I  x  I  -*'*>  I  x  I  -*® 

integrated  twice  by  parts.  The  second  term  is  the  Fourier  transform 
of  R(x,t),  i.e.,  R(A,t).  We  get  the  following: 


w(  A , t )  =  -  kA  w(A,t)  +  R( A , t) 


w( A , 0)  =  f ( A )  , 


f(A)  =  e  f(x)  dx 


This  is  an  ordinary  differential  equation  for  w(A,t).  Its  solu¬ 


tion  is 


,2,  f  .2,  *  .  2,  - 

..  -A  kt  A  ks  .  ,  -A  kt  ..... 

w(A,t)  =  e  e  R(A,s)  ds  +  e  f(A)  . 


Consequently,  inverting  the  Fourier  transform,  we  have: 
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t 


»(x,c)  -  JJ 


r°° 

o 

iAx 

e 

-A  kt 
e 

J-oo 

pOO 

>  -  h 

iAx 

e  e 

-00 

.  2.  .  2.  _  • 

A  ks  -A  kt  .... 

e  R(A,s)ds+e  f(A) 


dA 


w(x,  t)  =  7T- 


f(A)  dA  +  ^ 


iAx-A  kt 


.  2,  A 

e  R(A,s)  ds  dA 


The  convolution  theorem  for  Fourier  transforms  gives: 


-QO 


2n 


*  \  .  A  A 

e1  X  g(A ; t)  h ( A ; t )  dA  *  g  *  h(x,t) 


r00 

g(4.t)  h(x  -  f,t)  d? 

-00 


and  therefore, 


2tt 


r00 

iAx 

- A^kt  2 , . , 

dA  = 

e 

e  f  (A) 

''-oo 

. 

df  f(0 


e-[  (x-F)2Akt] 

y^ffkt 


The  second  integral, 
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<3T  4  32T  1  Teo 

9t  t  2  r_  r_ 

T  ax  T  T 


E  (x.t) 


T(x , o)  =  T 


the  solution  of  which  is,  by  Eq .  (A. 2): 


T(x.t)  =  e 


■<t/rT) 


0  J  -00 


e-[(x-0  /4k(  t-s )  ]  rT  Teo  e2( 

-  - -  e  •  -  1  +  - r - 

-  s)  rT  E 


+  d? 


r[(x-A)V4ktl 

JUnut 


(A. 3) 


h  *4 

where  k  -  —  . 

rT 


Equation  (A. 3),  by  virtue  of  the  identity, 


-[(x-z)2/u2J 
dz  e  =  Jir  u 


is  equivalent  to 
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Appendix  B 

ASYMPTOTIC  LIMIT  TO  EQ.  (7) 


We  analyze  the  asymptotic  behavior  of  the  expression 
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/W  f. 


f00 

Q(0- 

liin 

u-*» 

<X2 

r  u  (x  - 

u  expr  7z  -  -4^J 


JkKXX 


du 


df  ,  (B.3) 


now 


r  U  (x  -  0^1 

r  expr  ^ '  -5sr- J 


74. 


du  = 


KU 


rU  0 

[  4xu  (x-£)  1  , 

eXp[-  kicrZ  ~  ~ 4xu  J  dU  ‘ 


We  let  u'  =  4*eu  in  Eq.  (B.4)  and  obtain 


( B  .4) 


r  u' 

L  r  l  ^  ~  -=H 

“  Jo  JZ7 


du'  . 


(B.5) 


Now  set  u'  =  w  ,  du'  =  2w  du  =  27u'  dw.  Then  Eq.  (B.5)  becomes: 


[v/4  itu 


JL^ 

4»c 


exp 


2  2 
w  o «-er 


4  «r  2 

T  w 


dw 


(B.6) 


From  Abramowitz  and  Stegun  [15],  we  know  that: 


I 


.2  2  .,  2  2.  . 

e(-a  w  ‘  <b  /w  )]  dw  _  7a 

4x 


2ab 

e  erf 


(aW  +  w  )  +  e’2ab  erf  (aW  •  l}\  + 


where 
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and , 


Hence,  from 


erf (z) 


a  = 


y  4 


KTn 


b  -  Vo2  -  ' 


iq .  (B.6)  we  get  the 


following : 


p74/cu 

2 

w 

(x-f)2 

exp 

•  0 

4/crT 

w2 

dw 


74 


*KT„ 


|x-$|/7«t 


erf 


.74 


_u _  |  x-£ 

■ -  w 


KT. 


|x-*|/7^ 


+  e 


74 


XU 


erf 

W  _  |x-$| 

.  74»crT  W 

}j 

0 

J*KT^ 

|x-||/7«rT 

erf 

u  +  Ml 

4 

e 

~r 

,/r^  7 4xu  , 

|x-f |/7«r 


+  e 


kzil 


74ku  J 


x-$|/7«t 


-  e 


erf  (+<»)  -  e 


|x-(f  \/Jkt 


erf  (-oo) 
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Jk*KT^ 

r. 

|x-$|/y*rT 

1  +  erf 

u 

ix-«r 

4 

r 

.J~T 

JUku  j 

.  erfc 

*  ' 

u  |x-P| 

' 

>T  y4«u  . 

"  » 

where  erfc(u)  =  1  -  erf(u) . 

As  we  take  lim,  u  -*  ®  and  the  above  expression  tends 
t-*0° 


Ju-KKT^  -(  |x-£  \/Jkt^) 

- - -  •  2  e 


Hence  the  expression  we  seek  is: 


lim 

U~*oo 


r 

Jo 


exp 


u  (x  -  Q 


2-, 


r_  4*cu 
T 


JUku 


dw 


2_  ..  fA"u  f.  _»L  .  (X  -  <)2 

JQ  P[  4«rT  v2 


2  JUkkt^  "(|x-£| /Jkt^) 

Tk  •  -4 - 2  e 


JkKKT^  -(  |x-£|/ArT) 
4x  6 


We  substitute  Eq.  (B.9)  back  into  Eq .  (B.3)  to  get: 


|x-f \/Jkt^ 


(B.7) 


toward 


(B.8) 


dw 


(B.9) 
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- -  —CO 

-  i 

/n  t  _/c 
'  T  oo 


Q(0  e 


-( i*-£  i/y*o 


i  , 

/'"l  J-co 


Q(0  e 


-(|x-f \/Jkt  ) 


Hence,  we  find  that 


lim  T(x , t)  =  T  1  +  — 

t-®  C°l  2/ 


r°o 

1 _ 

X  J-00 


QCO  e 


( |x-f \/Jkt  ) 


It  is  then  clear,  as  we  consider  the  limiting  behavior  as  k  -*  0  (i.e., 
no  diffusion) ,  that 


lim  T(x, t) 

t-*00 
K-*  0 


T  ♦  ^ 

60  Ep2 


since  the  kernel 


-( |x-£ \/Jkt  ) 


becomes  a  6-function  as  k  -*  0. 
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Appendix  C 

STEADY- STATE  SOLUTION  TO  EQ.  (18) 


The  steady-state  form  of  Eq .  (18)  is: 


2  d2n  2  N0  dATe  N0  AT 

Hi  2  +  RTLN  T  2  n  7  T  ATe  ' 

dx  eo  dx  eo 


for  boundary  conditions  lim  n  =  0.  This  is  of  the  form 

x-»±oo 


y"  -  a  y  =  - f (x)  y  (±°°)  =  0  , 


The  solution  is  found  by  the  method  of  variation  of  parameters  by 
assuming  y  =  meax  +  ^2e_ax-  We  solve  the  system  thus: 

ax 

yl  =  6  ’ 

-ax 

Yo  =  e 


"I  yi  +  *2  y2 


M1  yl  +  U2  y2 


=  -f(x)  , 


and  find 


,  1  axc.  , 

m2  2a  6  f 


,  1  -ax,,  . 

"l  “  -  21  e  f(x) 


whence 
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ax 

e 

2a 


rx 


-as  t/  .  .  .  e 

e  f(s)  ds  + 


-ax 


2a 


eaS  f(s)  ds , 


Putting  £1  =  +<*>  and  ^2  =  allows  lim  y  =  0,  and  hence, 


Y  =  oZ 


2a 


-a(x-s)t/  .  ,  1 

e  f(s)  ds  +  2a 


-a(s-x) c .  .  , 

e  f(s)  ds 


r® 


Jl. 

2a 


-alx-s I 
e  1  1 


f(s)  ds  . 


^■00 


Now, 


-f(x)  = 


44  T  .2 
eo  dx 


*  N°  ^  ^  AT 


eo 


Using 


2 

d  AT 


dx 


1 

*4 


I  AT  -  T  q|  , 
^  e  eo  j 


where  Q  =  E^ /E^  , 


we  have : 


No 

r  *4 

f(x)  - 

T 

eo 

7  +  Rt  — 

4 

ATe  RT  N0  2 

4 

4 
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We  calculate 


n 

N, 


1 

2ln 


'I'(x')  e 


[x-x’ 

S 


dx' 


where  <t(x'  )  = 


eo 


7  +  R 


T  2 

S 


S 

ATe  "  RT  “2  Q 

4 


AT 

and  (x'  )  =  jj— 

eo  T 


U'-vl 

Q(y)  e  LT  dy 


(i.e.,  the  steady-state  temperature  solution  from  Appendix  B) .  Ex 
panding,  we  have 


't(x'  )  = 


S  1 

7  ^  ^ 


.  JjslzxL 

Q(y)  e  S  dy  - 


Ll 

rt  72  Q(*'> 


S’ 


7 

2S 


Q(y)  e 


_  l*'-yl 

S 


.  vi  - 

dy  +  — r— 

2S 


-  lx'-yl 

Q(y)  e  1^.  dy  -  RT 


7 

2S 


_  lx'-yl 

Q(y)  e  S  dy  + 


i 

2s 


Q(y)  e 


U'-yl 
S  dy 


then 


Q(x') 
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i-fhol^r 


A  lengthy  but  elementary  calculation  shows  that 


Thus 
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n_ 

N, 


2lnlt 


dy  Q(y ) 


ln  e 


Liizvj 

h 


*izl 


-  LT  e  T 


SS  2lnlt 


3  2  2 

4Lt  ln  -  lt 


dy  Q(y) 


x-y 

S  6  ^ 


L*=zl 


Lt  e  T 


2Lt2 


x-y 

dy  Q(y)  e  S 


0  2 


7Vt 

pOO 

dy  Q(y) 

-00 

e(-|x-y|/LN) 

(-|x-y|/LT) 

2  2 
s  -  s  J 

lt 

ln  J 

.  . 

^  e(-|x-y|/LN)  _  ^  e(-|x-y|/LT) 


V* 


A 


¥| 

2I4 


The  second  term  is  equal  to: 


dy  Q(y) 


w 

2L? 


2  2 

VS 


SS 


24 


l*-yl  r  t  2 

h  -  w 

2lt 


2  2 

VS 


W 


HA]  J 


dy  Q(y) 


S.  [sK]  rtH. 


2  , 2 


S'  SS  S 


_  Jbiiii 

■  h,  -i- 


And  so  we  find  the  following: 


dy  Q(y) 

|x-y|/LN) 

.  l*-yl 

e  S 

1 

s  . 
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4  4~4  44  4 

2  2  2  2 

4  44  4  4 


Hence , 


appears  to  be  singular  when  Lf|  =  Lj,  which,  according  to  Gurevich  [1]. 
occurs  around  300  km.  In  fact,  the  expression  has  a  removable 
singularity  here  as  we  shall  now  show.  We  let  L  =  and  10  =  Lj, 

then : 
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Lim 

1-1 


lira 

1-1 


1 


”  *-y 

e  1 


In  the  above  G(x,y,l,l0)  =  f(l)/g(l),  where  both  f(l)-0  and  g(l)-0 

1-1„.  Hence,  by  L'Hospital's  rule,  lim  G(l)  =  lim  f'(l)/g'(l) 

l^lo  l~lo 


f(l) 


-  [x-y[ 

e  1  - 

K 

V’ 

i 

1 

0 

lx-z_ 


f'(l) 


e 


x-y  | 
1 


+ 


i 


g'  (i) 


Therefore,  we  find  that 


1  im 
1-1 


f'(D 

g'(l) 


1 

4 


and , 


-  | x-y  r 
1 

o 


as 

,  where 
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li»  - 

S"'lt 


and  where  =  L^,  (at  300  km)  , 


which  is  well-defined. 
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Appendix  D 

SOLUTION  FOR  SQUARE  HEATED  REGION 


We  consider  an  idealized  joule-heating  source, 


E(x,t)  = 


E  -  L  <  x  <  L  and  o  <  t  <  t, 

o  1 

0  elsewhere  , 


and  integrate  the  expression 


^  i  r  d=  r 


where  k  =  L^/rT  • 


Thus  we  find: 


(t  -  s)  _  (x  -  Q 
A/c ( t  -s) 

Jhnic  ( t  -  s ) 


^  (f,s)  ds 
E 

P 


AT  ( x  ,  t )  _L_ 

T  tr2  rT 

eo  E  T 

P 


exp 

(t  -  s) 

rL  r 

rT 

.  _  L  <*  - 

74™  (t  -  s) 


4/c(t  -  s) 


d/:  ds 


where 


t  if  t  <  t^  , 
t^  if  t  >  t^  . 


We  shall  see  that,  in  the  event  the  source  is  turned  off  at  a  finite 
time  (t^  <  +o°) ,  the  perturbation  will  naturally  decay  to  zero.  Now. 
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rL 

0 

>L-x 

■ 

exp 

•-L 

(x  -  r  > 

d*  - 

exp 

-<L+x) 

Hi 

4<c  ( t  -  s ) 

4*(t  -  s) 

• 

-  ^ 

»  ■> 

' 

JknK ( t  -  s ) 

erf 

L  -  x 

+  erf 

L  +  x 

2 

/4x(t  -  s ) A 

/4x(t  -  s ) j 

where 


erf  (Z) 


dt  . 


Hence,  we  find: 


AT 
_ e 

T 

eo 


.2  2 


rT 

exp 

•  o 


*  ' 

L  -  x 

+  erf 

■ 

L  +  x 

./4x(t  -  s) 

./4x(t  -  s)J 

ds  . 


We  then  use  the  following,  from  Abramowitz  and  Stegun  [15]: 


j"  eav  erf  /b/v  dv  =  ^  jeaV  erf  Jb/v  +  i 


[ av- (b/v) ] 


|w ( 7 av 


+  i  Jb/v)  +  w(-  J av  +  i 


/b/v)]} 


2 

where  w(z)  -  e  erfc(-iz), 

erfc(z)  ■  1  -  erf (z) , 
i  =■  /- 1) . 
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Hence,  we  have  the  following: 


rt-T 


exp 


(t  -  s) 


erf 


L  ±  x 


1/4*  (t  -  s) 


ds  =  - 


exp 


-u 


erf 


[L  ±  x 


J  4 


XU 


du 


fT  ' 

[  exp 

~U] 

erf 

L  ±  x 

I 

l 

Tt\ 

.  y4*(t  -  s>_ 

+  2  exP 


u_  _  (L  ±  x)‘ 
t„  4xu 


y-u/r. 


.  L  ±  x 

+  l  J - L 


JUk  u  J 


+  w 


-  y-u/t  +  i 


L  ±  x 

1  J - ! - L 


y4/cu  J 


u=t-T 


u=t 


Note,  since  u  >  o,  ±J-u/t ^  =  ±i  Ju/t^  . 


Now  we  find  that 


therefore , 
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If  t  <  t  T  =  t,  and  observing  that  lim  erfc(u)  =  o  and  lim  erf(u)  = 

1  li-«»  U-*ao 

1,  we  find  that 
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If  the  heat  source  is  never  turned  off,  that  is,  taking  t^  =  +  ®, 
we  have : 


This  is  in  agreement,  of  course,  with  the  steady-state  solution  ob¬ 
tained  by  using  the  formula  in  Appendix  B: 


where 

|x|  <  L  . 

!xl  >  L  . 

Moreover,  if  ®  and  we  take  T  =  t]^  in  Eq .  (D.l),  the  limit  is 

clearly  zero. 
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